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Abstract 

The Relativistic Hartree Bogoliubov theory is applied in the mean- 
field approximation to the description of properties of light nuclei with 
large neutron excess. Pairing correlations and the coupling to particle 
continuum states are described by finite range two-body forces. Using 
standard parameter sets for the mean-field Lagrangian and the pairing 
interaction of Gogny-type, self-consistent solutions in coordinate space 
are calculated for the ground states of a number of neutron-rich nuclei. 
The model predicts the location of the neutron drip-line, reduction of 
the spin-orbit interaction, rms radii, changes in surface properties, 
the formation of neutron skin and of neutron halo. 
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1 Introduction 



Experimental and theoretical studies of exotic nuclei with extreme isospin 
values present one of the most active areas of research in nuclear physics. Ex- 
periments with radioactive nuclear beams provide the opportunity to study 
very short-lived nuclei with extreme neutron to proton ratio N/Z. New theo- 
retical models and techniques are being developed in order to describe unique 
phenomena in nuclei very different from those usually encountered along the 
line of stability. On the neutron-rich side in particular, exotic phenomena 
include the weak binding of the outermost neutrons, pronounced effects of 
coupling between bound states and the particle continuum, regions of neu- 
tron halos with very diffuse neutron densities and major modifications if the 
shell structures. Neutron drip lines of relatively light nuclei have become ac- 
cessible in experiments, and properties of these exotic objects are currently 
being studied with a variety of theoretical approaches. Because of their rel- 
evance to the r-process in nucleosynthesis, nuclei close to the neutron drip 
line are also very important in nuclear astrophysics. Detailed knowledge of 
their structure and properties would help the determination of astrophysical 
conditions for the formation of neutron rich stable isotopes. 

The main problem in the theoretical description of drip-line nuclei arises 
from the closeness of the Fermi level to the particle continuum: particle- 
hole and pair excitations reach the continuum. The coupling between bound 
states and the particle continuum has to be explicitly included in any model 
description. The Relativistic Hartree Bogoliubov (RHB) theory, which is a 
relativistic extension of the Hartree Fock Bogoliubov theory, provides a uni- 
fied description of mean-field and pairing correlations. A fully self-consistent 
RHB theory in coordinate space correctly describes the coupling between 
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bound and continuum states. The theory provides a framework for describ- 
ing the nuclear many-body problem as a relativistic system of baryons and 
mesons. In the self-consistent mean-field approximation, the theory has been 
used to perform detailed calculations for a variety of nuclear structure phe- 
nomena, as well as for the dynamics of heavy-ion collisions. Although con- 
ventional non-relativistic Hartree-Fock models and the relativistic mean-field 
theory predict very similar properties for nuclei close to the /3-stability line, 
recent applications to the structure of drip-line nuclei have shown significant 
differences, especially in the isospin dependence of the spin-orbit potential. 

In Refs. [[I], |H we have investigated, in the framework of relativistic 
Hartree-Bogoliubov theory, light nuclear systems with large neutron excess. 
The isospin dependence of the spin-orbit interaction in light neutron-rich 
nuclei has been studied in Ref. [0]. In |] we have described the micro- 
scopic details of the formation of neutron halo in the mass region above the 
s — d shell. Pairing correlations and the coupling to particle continuum states 
have been described by finite range two-body Gogny-type interaction. Fi- 
nite element methods have been used in the coordinate space discretization 
of the coupled system of Dirac-Hartree-Bogoliubov and Klein-Gordon equa- 
tions P, §]. Calculations have been performed for the isotopic chains of Ne 
and C nuclei. We find evidence for the occurrence of multi-neutron halo in 
heavier Ne isotopes. The detailed microscopic analysis has shown how the 
properties of the lf-2p orbitals near the Fermi level and the neutron pairing 
interaction determine the formation of the halo. On the other hand our cal- 
culations have shown that, although for Ne the neutron drip was found at 
N = 32, for C already the N = 20 isotope is unbound. Preliminary results 
for Oxygen have placed the neutron drip at iV = 18 or N = 20, depending 
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on the effective force used in the mean-field calculations. This is a very in- 
teresting result: by adding just one or two protons, the resulting mean- fields 
bind up to twelve more neutrons. 

In the present work we extend our investigation of Ref . |]TJ , and study in 
detail the neutron drip-line for the elements C, N, O, F, Ne, Na and Mg. In 
the region around N = 20 we expect that these nuclei are to a good approx- 
imation spherical, or that an eventual small deformation would not change 
our results quantitatively. In particular, we study how the proton number af- 
fects the neutron single-particle levels around N = 20 and the position of the 
neutron drip. We also discuss changes in the spin-orbit interaction, surface 
difluseness and the formation of the neutron halo. In Sec. 2 we present an 
outline of the relativistic Hartree-Bogoliubov theory and discuss the various 
approximations. In Sec. 3 the theory is applied to the calculation of ground 
state properties of a number of light neutron-rich nuclei. A summary of our 
results is presented in Sec. 4. 

2 The relativistic Hartree-Bogoliubov model 

Detailed properties of nuclear matter and finite nuclei along the /3-stability 
line have been very successfully described in the framework of relativis- 
tic mean field models 0, || |7|, ||. In comparison with conventional non- 
relativistic approaches, relativistic models explicitly include mesonic degrees 
of freedom and describe the nucleons as Dirac particles. Nucleons interact in 
a relativistic covariant manner through the exchange of virtual mesons: the 
isoscalar scalar cr-meson, the isoscalar vector a»-meson and the isovector vec- 
tor p-meson. The model is based on the one boson exchange description 
of the nucleon-nucleon interaction. We start from the effective Lagrangian 
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density 

£ = if) («7 ■ d — m) if) 

+ \{daf - U(a) - ifl^ + \ml^ - + \m\f - -Y ^ 

-g^ifxjif) - g^l ■ ojif) - g p if)^ ■ pfip - e^7 • A ^ 1 n - if) . (1) 

Vectors in isospin space are denoted by arrows, and bold-faced symbols will 
indicate vectors in ordinary three-dimensional space. The Dirac spinor if) 
denotes the nucleon with mass m. m a , m u , and m p are the masses of the 
cr-meson, the cu-meson, and the p-meson. g a , g u , and g p are the correspond- 
ing coupling constants for the mesons to the nucleon. e 2 /47r = 1/137.036. 
Coupling constants and unknown meson masses are parameters, adjusted to 
fit data on nuclear matter and finite nuclei. U(a) denotes the non-linear a 
self-interaction || 

U(a) = l -mW + \<hP* + \g^\ (2) 

and f^", R tJ,v , and F^ u are field tensors 

= (3) 
RT = dPf-Vp* (4) 
= d»A v -d v A". (5) 

The lowest order of the quantum field theory is the mean-field approximation: 
the meson field operators are replaced by their expectation values. The A 
nucleons, described by a Slater determinant |<&) of single-particle spinors 
ipi, (i = 1,2, A), move independently in the classical meson fields. The 
sources of the meson fields are defined by the nucleon densities and currents. 
The ground state of a nucleus is described by the stationary self-consistent 
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solution of the coupled system of Dirac and Klein-Gordon equations. In the 
static case for an even-even system, there will be no currents in the nucleus, 
and the spatial vector components u), p 3 and A vanish. The Dirac equation 
reads 

| -ia ■ V + P(m + g a a) + g^u + g p T 3 p 3 + e^—^-^A° j fa = e^i (6) 
The effective mass m*(r) is defined 

m* (r) = m + g a cr(r), (7) 

and the potential V(r) 

V(v) = g uU °{r) +g p T 3 p 3 (v) + e^^-A°(v). (8) 

In order to describe ground-state properties of spherical open-shell nuclei, 
pairing correlations have to be taken into account. For nuclei close to the 
/3-stability line, pairing has been included in the relativistic mean-field model 



in the form of a simple BCS approximation QlO | . However, for nuclei far from 
stability the BCS model presents only a poor approximation. In particular, 
in drip-line nuclei the Fermi level is found close to the particle continuum. 
The lowest particle-hole or particle-particle modes are often embedded in the 
continuum, and the coupling between bound and continuum states has to be 
taken into account explicitly. The BCS model does not provide a correct 
description of the scattering of nucleonic pairs from bound states to the 
positive energy continuum O, [EJ. It leads to an unbound system, because 
levels in the continuum are partially occupied. Including the system in a box 
of finite size leads to unreliable predictions for nuclear radii depending on the 
size of this box. In the non-relativistic unified description of mean-field 

and pairing correlations is obtained in the framework of the Hartree-Fock- 



Bogoliubov (HFB) theory in coordinate space fil l. The ground state of a 



nucleus |$ > is represented as the vacuum with respect to independent quasi- 
particles. The quasi-particle operators are defined by a unitary Bogoliubov 
transformation of the single-nucleon creation and annihilation operators. The 
generalized single-particle Hamiltonian of HFB theory contains two average 
potentials: the self-consistent field f which encloses all the long range ph 
correlations, and a pairing field A which sums up the pp-correlations. The 
expectation value of the nuclear Hamiltonian < §\H |$ > can be expressed as 
a function of the hermitian density matrix p, and the antisymmetric pairing 
tensor k. The variation of the energy functional with respect to p and k 



produces the single quasi-particle Hartree-Fock-Bogoliubov equations [13 



h-X A \ fU k \ „ (U k 



A* -h + Xj\V k J- Ek {vJ- (9) 
HFB-theory, being a variational approximation, results in a violation of basic 
symmetries of the nuclear system, among which the most important is the 
non conservation of the number of particles. In order that the expectation 
value of the particle number operator in the ground state equals the number 
of nucleons, equations (0) contain a chemical potential A which has to be 
determined by the particle number subsidiary condition. The column vec- 
tors denote the quasi-particle wave functions, and E k are the quasi-particle 
energies. 

If the pairing field A is a constant, HFB reduces to the BCS-approximation. 
The lower and upper components U k (r) and V k (r) are proportional, with the 
BCS-occupation amplitudes as proportionality constants. For a more gen- 
eral pairing interaction this will no longer be the case. As opposed to the 
functions U k (r), the lower components V k (r) are localized functions of r, as 
long as the chemical potential A is below the continuum limit. Since the den- 
sities are bilinear products of V k {r) (see Eqs.([T5|)- ([T8|) ), the system is always 
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localized. According to the theorem of Bloch and Messiah |R|, an HFB wave 
function can be written either in the quasiparticle basis as a product of inde- 
pendent quasi-particle states, or in the canonical basis as a highly correlated 
BCS-state. In the canonical basis nucleons occupy single-particle states. If 
the chemical potential is close to the continuum, the pairing interaction scat- 
ters pairs of nucleons into continuum states. Because of additional pairing 
correlations, particles which occupy those levels cannot evaporate. Mathe- 
matically this is expressed by an additional non-trivial transformation which 
connects the particle operators in the canonical basis to the quasi-particles 
of the HFB wave function. 



The relativistic extension of the HFB theory was introduced in Ref. [[14 
In the Hartree approximation for the self-consistent mean field, the Rela- 
tivistic Hartree-Bogoliubov (RHB) equations read 

fh D -m-X A \{U k (r)\_ (U k (v)\ 

{ -A* -h D + m + Xj \Vk(v)J -* k \V k (r)J- 

where ho is the single- nucleon Dirac Hamiltonian(^), and m is the nucleon 

mass. The RHB equations have to be solved self-consistently, with potentials 

determined in the mean-field approximation from solutions of Klein-Gordon 

equations 

[-A + ml\a(r) = - g<7 p s (r) - g 2 a 2 (r) - g 3 a 3 (r) (11) 

|-A+m>°(r) = ~g u Pv{v) (12) 
[-A + m>°(r) = -g pP3 (r) (13) 
- AA°(r) = e Pp (r). (14) 

for the sigma meson, omega meson, rho meson and photon field, respectively. 
The spatial components u, p, and A vanish due to time reversal symmetry. 
Because of charge conservation, only the 3-component of the isovector rho 



meson contributes. The source terms in equations flTT|) to ([14]) are sums of 
bilinear products of baryon amplitudes 

Ps(r) = E ^(r)A(r), (15) 

E k >0 

Pv(r) = E Vi(r)V k (r), (16) 

E k >0 

PaW = E v£(r) T3 V k (v), (17) 



E k >0 



p em (r) = E V2(r) l -^V k {v). (18) 

E k >0 Z 

where the sums run over all positive energy states. For M degrees of free- 
dom, for example number of nodes on a radial mesh, the HB equations are 
2M-dimensional and have 2M eigenvalues and eigenvectors. To each eigen- 
vector (Uk, Vjfc) with eigenvalue E k , there corresponds an eigenvector (V k *, C/|) 
with eigenvalue — E k . Since baryon quasi-particle operators satisfy fermion 
commutation relations, it is forbidden to occupy the levels E k and — E^ si- 
multaneously. Usually one chooses the M positive eigenvalues E k for the 
solution that corresponds to a ground state of a nucleus with even particle 
number. 

The system of equations fllPD , and dTTf) to (|HD, is solved self-consistently 



in coordinate space by discretization on the finite element mesh || ^| . In the 



coordinate space representation of the pairing field A in (10), the kernel of 
the integral operator is 

A a6 (r, r') = ~E Vabcdir, r> C(i (r, r'). (19) 

c,d 

where a, b, c, d denote all quantum numbers, apart from the coordinate r, 
that specify the single-nucleon states. V a b c d{r, r') are matrix elements of a 
general two-body pairing interaction, and the pairing tensor is defined as 

K cd (r,r')= E U: k {r)V dk {v'). (20) 

E k >0 

9 



The integral operator A acts on the wave function Vfc(r): 



(A^)(r) = J] / <fVA a6 (r,r')^(r') 



(21) 
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The eigensolutions of Eq. ([TOj) form a set of orthogonal (normalized) sin- 
gle quasi-particle states. The corresponding eigenvalues are the single quasi- 
particle energies. In the self-consistent iteration procedure we work in the ba- 
sis of quasi-particle states. The self-consistent quasi-particle eigenspectrum is 
then transformed into the canonical basis of single-particle states. The canon- 
ical basis is defined to be the one in which the matrix R^i = (Vfc(r)|Vfc/(r)) is 
diagonal. The transformation to the canonical basis determines the energies 
and occupation probabilities of single-particle states, that correspond to the 
self-consistent solution for the ground state of a nucleus. 

For nuclear systems with spherical symmetry the fields cr(r), u>°(r), p°(r), 
and A°(r) depend only on the radial coordinate r. The nucleon spinors Uk 
(V)fc) in (|T0D are characterized by the angular momentum j, its ^-projection 
m, parity tc and the isospin £ 3 = ±| for neutron and proton. The two Dirac 
spinors C/&(r) and 14 (r) form a super-spinor 



g(r) and f(r) are radial amplitudes, \ T is the isospin function, the orbital 
angular momenta I and I are determined by j and the parity n 




(22) 



where 




(23) 




for 
for 




(24) 



and 




for 
for 




(25) 
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Vtji m is the tensor product of the orbital and spin functions 

£lj,i, m (0,(p,s)= ^^slrni\jm)xk m8 Y im{ & ^)- (26) 

m s ,mi 

It is useful to define a single angular quantum number k as the eigenvalue of 
the operator (1 + a ■ 1) 

(1 + a ■ 1) Q K>m = -Ktt K>m , (27) 

k = ±(j + 1/2) for j = I =f 1/2. (28) 

k = ±1, ±2, ±3, .... The equations for the radial amplitudes gu(v)( r ) an d 
fu(v){ r ) are derived from Eq. flTUp. The radial single-particle Dirac Hamil- 



tonian reads 

h D (r) = <x 3 ai(«9 r + r- 1 ) - a^r' 1 + a 3 m + a 3 S(r) + V (r), (29) 

where the scalar and vector potentials are 

S(r)=g a a(r), and V°(r) = g u u°{r) +g p r 3 p° 3 (r) + e^^-A°{r). (30) 

o"j (i=l,2,3) are the Pauli matrices. The integral operator of the pairing 
interaction takes the form 

CO 

A(r) = J dr'r' 2 A(r,r') (31) 
o 

The kernel of the integral operator is defined 

A^(r,r') = \Y:{™SAV\t%^jm«**M ( 32 ) 

a,a' 

where r and r' denote radial coordinates, a, a', a and a' are quantum num- 
bers that completely specify single- nucleon states: (n,l,j,m) or (n, K,m), J 
and M are the total angular momentum of the pair, and its z-projection, 
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respectively. In the particle-particle (pp) channel the pairing interaction is 
approximated by a two-body finite range Gogny interaction 

V pp (l,2) = e- ((ri_r2)/w)2 (Wi + BiP' 7 - HiP T - MiP a P T ), (33) 

i=l,2 

with parameters fii, Wi, Bi, Hi and Mi (i = 1,2). For the J = contribution 
to the pairing matrix elements, the kernel of the integral operator reads 

A K (ry) = ^V K t°(r,r>)K~ K (ry). (34) 

The pairing tensor is calculated 




where, for a quantum number k, the sum runs over all solutions in the spec- 
ified energy interval < E < E mSLX . If we define the Dirac spinors 




the radial Dirac-Hartree-Bogoliubov equations read 

(h D {r) -m- \)<$>u(r) + / dr'r' 2 A(r, r')<f> v (r') = E^ v (r) 

Jo 

(-h D (r) + m + A)$ y (r) + / dr'r' 2 A(r, r')$u(r') = E$ v (r) (37) 

Jo 

The meson and photon fields are solution of the Klein-Gordon equations 

{-d 2 r --d r + ml)a{r) = -g c p s (r) - g 2 a 2 (r) - g 3 a 3 (r) (38) 
r 

{-d 2 r -- r d r + ml)u\r) = -g uPv (r) (39) 

(-d 2 -\ + m 2 p )p°(r) = -g pP3 (r) (40) 

(-d 2 --d r )A°(r) = ep p (r). (41) 
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3 Neutron drip-line in light nuclei 

The details of the neutron drip-line will depend on the single-particle levels, 
position of the Fermi energy, pairing interaction and coupling between bound 
and continuum states. In order to ascertain that our microscopic predictions 
do not crucially depend on the choice of the effective interaction, we perform 
calculations for three parameter sets of the mean-field Lagrangian: NL1 [|15|1 , 
NL3 p|, and NL-SH |TJ. The effective forces NL1 and NL-SH have been 



frequently used to calculate properties of nuclear matter and of finite nuclei, 
and have become standard parameterizations for relativistic mean-field cal- 
culations. The parameter set NL3 has been derived recently [16j by fitting 
ground state properties of a large number of spherical nuclei. Properties 
calculated with the NL3 effective interaction are found to be in very good 
agreement with experimental data for nuclei at and away from the line of 
/3-st ability. 

In the pairing channel one should in principle use a one-meson exchange 
interaction V a bcd derived by the elimination of the mesonic degrees of freedom 
in the Lagrangian ([]]). However, as it was shown in Ref. [|14|], the standard 



parameter sets of the mean-field approximation lead to completely unrealistic 
pairing matrix elements. These parameters do not reproduce the scattering 
data. Since at present a relativistic pairing force with parameters consistent 
with those of the mean-field potential is not available, we use a phenomeno- 
logical Gogny-type finite range interaction in the pp-channel, a procedure 
which requires no cut-off and which provides a very reliable description of 
pairing properties in finite nuclei. The parameter set D1S [18| is used for the 
finite range Gogny type interaction. 

In Ref. we have studied in detail the formation of the neutron halo 

13 



in Ne isotopes. In order to be more complete and as an introduction to the 
discussion that follows, we briefly review the main arguments and results of 
that work. In Fig. la the rms radii for Ne isotopes are plotted as functions 
of neutron number. We display neutron and proton rms radii, and the N 1//3 
curve normalized so that it coincides with the neutron radius in 20 Ne. The 
neutron rms radii follow the mean-field N 1 / 3 curve up to N « 22. For larger 
values of N the neutron radii display a sharp increase, while the proton radii 
stay practically constant. This sudden increase in neutron rms radii has 
been interpreted as evidence for the formation of a multi-particle halo. The 
phenomenon was also observed in the plot of proton and neutron density 
distributions. The proton density profiles do not change with the number of 
neutrons, while the neutron density distributions display an abrupt change 
between 30 Ne and 32 Ne. The microscopic origin of the neutron halo has been 
found in a delicate balance of the self-consistent mean-field and the pairing 
field. This is shown in Fig. lb, where we display the neutron single-particle 
states lfr/2, 2p 3 / 2 and 2pi/ 2 in the canonical basis, and the Fermi energy as 
function of the neutron number. For N< 22 the triplet of states is high in the 
continuum, and the Fermi level uniformly increases toward zero. The triplet 
approaches zero energy, and a gap is formed between these states and all 
other states in the continuum. The shell structure dramatically changes at 
N> 22. Between N = 22 and N = 32 the Fermi level is practically constant 
and very close to the continuum. The addition of neutrons in this region of 
the drip does not increase the binding. Only the spatial extension of neutron 
distribution displays an increase. The formation of the neutron halo is related 
to the quasi-degeneracy of the triplet of states lf7/2, 2p 3 / 2 and 2p 1( / 2 . The 
pairing interaction promotes neutrons from the H7/2 orbital to the 2p levels. 
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Since these levels are so close in energy, the total binding energy does not 
change significantly. Due to their small centrifugal barrier, the 2p3/2 and 
2pi/2 orbitals form the halo. The last bound isotope is 40 Ne. For N > 32 
(A>42) the neutron Fermi level becomes positive, and heavier isotopes are 
not bound any more. A similar mechanism has been studied in Ref. [H5|, for 
the experimentally observed halo in the nucleus 11 Li. There the formation 
of the halo is determined by the pair of neutron levels lpi/2 and 2si/2- A 



giant halo has been also predicted for Zirconium isotopes 1 20| . In that case 
the halo originates from the neutron orbitals 2f 7 / 2 , 3p 3 / 2 and 3pi/ 2 . 

In Figs. 2, 3 and 4 we display the neutron single-particle energies in 
canonical basis for C, N, O, F, Ne, Na and Mg, in the region around neutron 
number N=20. In Figs. 2 and 3 the left hand side columns display results 
calculated with the NL1 effective interaction, NL3 results are shown in the 
central panels, and on the right hand side levels that correspond to the NL- 
SH parameterization are shown, as function of neutron number. For Carbon, 
in Fig. 4 we display results obtained with the NL3 and NL-SH parameter 
sets. The NL1 force does not produce a self-consistent solution for the ground 
states of C isotopes. We display the energy levels of the triplet of neutron 
states If 7/2, 2p3/2 and 2pi/ 2 , the neutron Fermi level and, as reference, the 
occupied orbital ld 3 / 2 . Other states in the continuum, as well as deep hole 
states, are too far away from the Fermi level and therefore not important 
in the present consideration. Results calculated with different effective in- 
teractions are rather similar, and in agreement with our conclusions for Ne 
isotopes: the location of the neutron drip is essentially determined by the 
triplet H7/2, 2p 3 / 2 and 2pi/ 2 . It appears that NL1 leads to strongest binding, 
followed by NL3 and NL-SH. For example, 28 is bound for NL1 and NL3, 
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but unbound for the NL-SH parameter set. As the neutron number N=20 is 
approached, for C, N, O and F the triplet is still high in the continuum, and 
the pairing interaction does not have enough strength to promote neutrons 
into these levels. Depending on the effective force, the neutron drip is found 
at N=18 or N=20. Starting from Ne, as we have already shown in Fig. 1, 
the triplet lf7/2, 2p 3 / 2 and 2pi/ 2 is much lower and neutrons above N=20 
begin to fill the lf 7 / 2 orbital. For Ne and Na, 2p 3 / 2 and 2pi/ 2 are close to 
the lf7/ 2 , and therefore they are also populated forming the neutron skin, 
or eventually the neutron halo. For Mg the neutrons above the s — d shell 
predominantly populate the lf 7 / 2 orbital, and one expects that deformation 
plays an important role in the precise location of the neutron drip. 

An interesting result is observed if one plots the neutron single-particle 
levels as function of the number of protons. This is done in Fig. 5 for 
N=18, N=20 and N=22 isotones. The levels 2p 3/2 , 2p 1/2 , lf 7/2 , and lf 5/2 , as 
well as the neutron Fermi level, are displayed for the three effective forces. 
For a given number of neutrons, increasing the proton number means going 
toward the /3-stability, i.e. away from the neutron drip-line. The binding 
increases, as reflected in the negative slope of the Fermi level. What we 
find very interesting is that the line that corresponds to the Fermi level is 
almost parallel to the lf 7 / 2 level. It appears that the binding of neutrons in 
these systems is determined by the position of the lf 7 / 2 orbital. The levels 
2p 3 / 2 , 2pi/2, and H5/2 change very little with Z. This is also true for lg9/ 2 , 
not shown in the figure, at an energy of ~ 20 MeV. On the other hand, the 
energy of the lf 7 / 2 orbital decreases considerably by adding more protons. If 
the spin-orbit splittings are compared, we notice that the energy difference 
for the doublet 2p 3 / 2 - 2pi/ 2 is almost constant as function of Z, while a very 
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strong isospin dependence is observed for lf 7 / 2 - H5/2. This is the crucial 
point and requires a careful explanation. In Ref. || we have analyzed the 
isospin dependence of the spin-orbit interaction for Ne and Mg isotopes. We 
have shown that the corresponding term in the potential is strongly reduced 
when going toward the drip-line. This reduction is especially pronounced 
in the surface region, and consequently one should observe a more diffuse 
surface. Since the densities of the / orbitals are located on the surface, while 
those of the p orbitals are contained in the bulk of the nucleus, the effect of 
reduction of the spin-orbit interaction will be much stronger for the lf 7 / 2 - 
If 5/2 doublet. 

In Fig. 6 we display the neutron densities for the N=20 isotones, cal- 
culated with the NL3 effective force. The figure corresponds to the central 
panel in Fig. 5. In the insert of Fig. 6 we include the corresponding values 
for the surface thickness and diffuseness parameter. The surface thickness 
s is defined to be the change in radius required to reduce p(r)/po from 0.9 
to 0.1 (po is the density in the center of the nucleus). The diffuseness pa- 
rameter a is determined by fitting the neutron density profiles to the Fermi 
distribution 

p(r) = Po (l + expC-^)) \ (42) 

where Rq is the half-density radius. In going away from the drip line, from 
28 to 32 Mg, the surface thickness decreases from 2.8 fm to 2.33 fm, and the 
diffuseness parameter decreases from 0.64 fm to 0.55 fm. We have shown 
the vector densities. However, the scalar densities display the same radial 
dependence. This means that the scalar potential S and the vector potential 
V behave in the same way in the surface region, i.e. both potentials display 
the same diffuseness. 
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These large changes in surface properties reflect the reduction of the spin- 
orbit term of the potential. In relativistic mean-field approximation, the spin- 
orbit potential originates from the addition of two large fields: the field of the 
vector mesons (short range repulsion), and the scalar field of the sigma meson 
(intermediate attraction). In the first order approximation, and assuming 
spherical symmetry, the spin orbit term can be written as 

K.o. = ~V h (r), (43) 
r or 



where Vi s is the spin-orbit potential M, |21 



V U = —{V-S). (44) 

m eff 

V and S denote the repulsive vector and the attractive scalar potentials, 
respectively. m e // is the effective mass 

m eff = m-l(V-S). (45) 

Using the vector and scalar potentials from the NL3 self-consistent ground- 
state solutions, we have computed from ( f43|) - fl45|) the spin-orbit terms for 
the N=20 isotones. They are displayed in Fig. 7 as function of the radial 
distance from the center of the nucleus. The magnitude of the spin-orbit 
term V s _ Qi increases as we add more protons, i.e. as we move away from 
the drip-line. In Fig. 8 we also plot the proton and neutron rms radii for 
the N=20 isotones, calculated with all three parameter sets. The results are 
very similar, and show how the neutron radii decrease in going toward the 
stability line. At the same time, of course, the proton radii increase. 

Finally, in Fig. 9 we compare the neutron rms radii for Na isotopes 



with recent experimental data from Ref . [22] . Theoretical values correspond 



to ground-state densities calculated with the NL3 effective interaction, and 



the D1S parameters for the Gogny interaction in the pairing channel. The 
blocking procedure has been used both for protons and neutrons. With the 
possible exceptions of N=ll, the theoretical values nicely reproduce the ex- 
perimental data, much better than the RMF parameter set that has been 



like in the case of Ne, no sudden change in the radii is observed. Instead 
of the neutron halo, what we observe in Na isotopes is more probably the 
formation of the neutron skin. In fact, in the insert of Fig. 9 we display the 
calculated and experimental differences between neutron and proton radii. 
The smooth increase of r n — r p after N=ll presents evidence for a large 
neutron skin. Of course, this does not prevent the formation of the halo 
for heavier Na isotopes. The discrepancy of the calculated and measured r n 
for N=ll probably comes from the fact that for this isotope the proton and 
neutron numbers are identical, and therefore both the odd proton and odd 
neutron occupy the s — d orbitals with the same probabilities. One should 
expect a short range interaction, possibly proton-neutron pairing. Such an 
interaction would lower the total energy, and the resulting radius should be 
smaller. An explicit proton-neutron interaction is not included in the theo- 
retical model, and therefore the calculated neutron radius is larger than the 
experimental value. 

4 Conclusions 

In the present work we have investigated properties of light nuclei with large 
neutron excess. Neutron-rich nuclei between Carbon and Magnesium have 
recently become available in experiments with radioactive nuclear beams, 
and therefore it is very important to study the properties that various the- 




The neutron rms radii display a uniform increase. Un- 
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oretical models predict for this region. We have calculated ground state 
properties using the Relativistic Hartree Bogoliubov theory. Based on the 
relativistic mean-field model on one side, and the HFB theory on the other, 
RHB provides a unified description of mean-field and pairing correlations. It 
describes the nucleus as a relativistic system of baryons and fermions. Pair- 
ing correlations are described by finite range two-body forces. Very little 
is known about the pp component of the effective interaction in relativistic 
models, and therefore we have used a finite range interaction of Gogny-type. 
Finite element methods have been used in the coordinate space discretiza- 
tion of the coupled system of Dirac-Hartree-Bogoliubov integro-differential 
eigenvalue equations, and Klein-Gordon equations for the meson fields. Co- 
ordinate space solutions are essential for a correct description of particle-hole 
and pair excitations into the continuum. For the mean-field Lagrangian we 
have used the well established effective interactions NL1, NL3 and NL-SH, 
and the D1S parameter set for the finite range Gogny type interaction in 
the pairing channel. In particular, we have investigated the location of the 
neutron drip-line as function of the number of protons. 

Model calculations have shown that the triplet of single-particle states 
near the neutron Fermi level: H7/2, 2p 3 / 2 and 2pi/ 2 , and the neutron pairing 
interaction determine the location of the neutron drip-line, the formation of 
the neutron skin, or eventually of the neutron halo. For C, N, O and F the 
triplet is still high in the continuum at N=20, and the pairing interaction is to 
weak to promote pairs of neutrons into these levels. All mean-field effective 
interactions predict similar results, and the neutron drip is found at N=18 or 
N=20. For Ne, Na, and Mg the states H7/2, 2p 3 / 2 and 2p 1 / 2 are much lower in 
energy, and for N> 20 the neutrons populate these levels. The neutron drip 
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can change by as much as twelve neutrons. The model predicts the formation 
of neutron skin, and eventually neutron halo in Ne and Na. This is due to the 
fact that the triplet of states is almost degenerate in energy for N> 20. For 
Mg the lf 7 / 2 lies deeper and neutrons above the s — d shell will exclusively 
populate this level, resulting in a deformation of the mean field. When we plot 
the neutron single-particle levels as function of the number of protons, it turns 
out that the binding of neutrons is determined by the H7/2 orbital which is 
found to be parallel to the slope of the Fermi level. The reduction of the spin- 
orbit splitting lf 7 / 2 - lf 5 / 2 close to the neutron-drip line has been related to the 
isospin dependence of the spin-orbit interaction. Density distributions have 
been analyzed, and the resulting surface thickness and diffuseness parameter 
reflect the reduction of the spin-orbit term of the effective potential in the 
surface region of neutron-rich nuclei. Model calculations reproduce the recent 
experimental data on the neutron rms radii for Na isotopes. 

The results of the present work clearly show that a correct theoretical 
description of nuclei with large neutron excess necessitates the following es- 
sential ingredients: 

i) self-consistent solution for the mean-field in order to obtain the correct 
increase of the surface diffuseness, 

ii) the correct isospin dependence of the spin-orbit term in the mean-field 
potential. At present only relativistic models provide reliable results for 
the characteristic spacings of single particle levels at the Fermi surface. 

iii) pairing correlations describe the scattering of nucleon pairs for states 
close to the Fermi surface. Pairing leads to the partial occupation of 
single particle levels with low orbital angular momentum. 
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iv) pairing correlations in the continuum should be described by RHB in 
coordinate space, with finite-range pairing interaction. In this proce- 
dure the scattering of nucleon pairs does not lead to unbound configu- 
rations. 

Thus Relativistic Hartree-Bogoliubov theory in coordinate space provides 
an appropriate framework in which properties of exotic nuclei on both sides of 
the valley of /3-stability can be studied. However, to calculate ground-state 
properties of heavier nuclei beyond Na or Mg, deformations of the mean- 
field have to be taken into account. The inclusion of deformation in the 
coordinate space formulation of RHB with finite range pairing still presents 
considerable difficulties. It is also clear that the details of the scattering 
processes around the Fermi surface, and the partial occupation of levels with 
low orbital angular momentum in the continuum, are determined by the 
strength and functional dependence of pairing correlations. As long as a fully 
microscopic derivation of an effective pairing in nuclei is not available, it is 
justified to use a phenomenological interaction, as for example the Gogny 
force. Another open problem is whether one needs a density dependent 
pairing interaction. 
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Figure Captions 

Fig.l Calculated proton and neutron rms radii for Ne isotopes (a), and the 
lf-2p single-particle neutron levels in the canonical basis (b). 

Fig. 2 The neutron single-particle levels ld 3 / 2 , H7/2, 2p3/2 and 2pi/ 2 for N, 
O, and F isotopes. Solid lines denote the neutron Fermi level. The 
energies correspond to ground-state solutions calculated with the NL1, 
NL3 and NL-SH effective forces of the mean-field Lagrangian. The 
parameter set D1S is used for the finite range Gogny-type interaction 
in the pairing channel. 

Fig. 3 Same as in Fig. 2, but for the isotopes of Ne, Na and Mg. 

Fig. 4 Same as in Fig. 2, but for C isotopes and the NL3 and NL-SH effective 
interactions. 

Fig. 5 The neutron single-particle levels 2p 3 / 2 , 2px/ 2 , lf 7 / 2 , and lf 5 / 2 for the 
N=18, N=20 and N=22 isotones, as function of the proton number. 
Solid lines denote the neutron Fermi level. The first, second and third 
row correspond to calculations performed with the NL1, NL3 and NL- 
SH effective forces, respectively. 

Fig. 6 Neutron density distributions for the N=20 isotones, calculated with 
the NL3 effective force. In the insert the corresponding values for the 
surface thickness and diffuseness parameter are included. 

Fig. 7 Radial dependence of the spin-orbit term in self-consistent solutions 
for the ground-states of the N=20 isotones, calculated with the NL3 
effective force. 
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Fig. 8 Proton and neutron rms radii for the N=20 isotones. The panels 
correspond to results obtained with the NL1, NL3 and NL-SH effective 
forces, respectively. 

Fig. 9 Comparison of the neutron rms radii for Na isotopes with experi- 
mental data from Ref. [35] . Theoretical values (squares) correspond 



to ground-state densities calculated with the NL3 effective interaction. 
The calculated differences between neutron and proton radii are dis- 
played in the insert. 
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